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ABSTRACT 

Using a cosmological N-body numerical simulation of the formation of a galaxy cluster- sized 
halo, we analyze the temporal evolution of its globular cluster population. We follow the dynam¬ 
ical evolution of 38 galactic dark matter halos orbiting in a galaxy cluster that at redshift z = 0 
has a virial mass of 1.71 x . In order to mimic both ’’blue” and ’’red” populations of 

globular clusters, for each galactic halo we select two different sets of particles at high redshift 
(z Ri 1), constrained by the condition that, at redshift z = 0, their average radial density profiles 
are similar to the observed profiles. As expected, the general galaxy cluster tidal field removes 
a significant fraction of the globular cluster populations to feed the intracluster population. On 
average, halos lost approximately 16% and 29% of their initial red and blue globular cluster 
populations, respectively. Our results suggest that these fractions strongly depend on the orbital 
trajectory of the galactic halo, specifically on the number of orbits and on the minimum pericen¬ 
tric distance to the galaxy cluster center that the halo has had. At a given time, these fractions 
also depend on the current clustercentric distance, just as observations show that the specific 
frequency of globular clusters Sn depends on their clustercentric distance. 

Subject headings: galaxies: clusters: general — galaxies: star clusters: general — methods: numerical 


1. Introduction 

Galaxy clusters are extreme environments 
where thousands of galaxies orbit inside a massive 
system of ^ a few Mpc size, and with 

a typical velocity dispersion of ^ 1000 km s“^. 
In these regions, galaxies evolve under the in¬ 
fluence of the strong general tidal field gener¬ 
ated by the clusters gravitational potential and 
also by individual mutual g ravitational interac¬ 
tions with other galaxies (e.g . iMerrittlflOSSl 1984 ; 
Malumuth fc Richstonel 19841 1. 


From the theoretical point of view, it is ex¬ 
pected that the global tidal field in a massive 
galaxy cluster will be strong enough to t runcate 
dark m atter halos of individual galaxies (Merritt 
Il983l l. IPiemand. Kuhlen fc MadaiJ (2007) found 
that the removal of the dark matter halo happens 
from the outside in, i.e. the outer particles of a 
halo are less bound and thus more prone to be re¬ 


moved. Although the baryonic material is much 
more concentrated than the dark matter, it can 
also suffer the effects of the galaxy cluster tidal 
field. Using N - body/ hydrodynamic simulations, 
Limousin et al. ( 2009ll investigated the effect of 
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tidal stripping on galaxy clusters, tracking both 
the dark matter halo and the stellar component, 
and found that the dark matter component is pref¬ 
erentially stripped, while the stellar component is 
less affected by tidal forces. 

Globular cluster systems are among the most 
extended galactic stellar components; they can 
reach out well beyond the detectable stellar halo 
and are also prone to being tidally stripped. Since 
dark matter density in central regions of galaxy 
clusters is expected to be higher than in their 
outskirts, this could result i n more compact ha - 
los in these central regions ( Bullock et al.ll^Oll) . 
A similar effect could be expected for globu¬ 
lar clusters, lowering their specific frequency Sn 
(number per luminosity unit) in the central re¬ 
gions o f galaxy clusters compared to th e out¬ 
skirts. iForbes. Brodie fc Grillm^ ( 1997 1 found 


1 







































evidence of this effect in four galaxies in the cen¬ 
tral region of the Fornax cluster with a marginal 
dependence of the spe cific frequency on the clus¬ 


tercentric distance. Coenda. Muriel fc DonzeUil 


(I 2 OO 9 II searched for similar evidence using data 
from the ACS Virgo Cluster Survey on board the 
Hubble Space Telescope. These authors analyzed 
the Sn in 13 elliptical galaxies and found that Sn 
increases as a function of both projected and 3D 
Virgo clustercentric distance, while the projected 
number density of background globular clusters 
decreases. These results are interpreted as di¬ 
rect evidence of tidal stripping of globular clusters 
due to the gravitational p otential of the galaxy 
cluste r. In a similar sense, lAlamo-Martmez et al 


(|2013ll studied the globular cluster system in the 
center of the massive galaxy cluster Abell 1689 
and estimated that 80, 000 out of a population 
of 162,000 are possibly part of the intrac luster 
medium (see also Webb. Sills fc Hai^ 2013ll . Ev¬ 
idence of tidal stripping of globular clusters has 
also bee n found in less den se environments ( see for 
instance iBlom et ^l2014h . Hudson et al.l (j2014h 
use galaxy halo masses based on weak lensing to 
compute the ratio between the total mass in glob¬ 
ular clusters and the halo mass. They found that 
this ratio is constant with an intrinsic scatter of 
0.2 dex and suggested that some of this scatter 
is due to differing degrees of tidal stripping of 
the globular cluster systems between central and 
satellite galaxies. 

In the last two decades or so, growing evidence 
has been collected of a bimodal color distribution 


of g l obular clusters (IGebhardt fc Kissler-Patia 


1999; Knndn fc Whitmorell200lHBrodie fc Stradei 


20061) . This is usually interpreted as the conse¬ 
quence of differences in metallicity between two 
co-existing globular cluster populations: one red 
and metal-rich [Fe/H] ^ —0.5 and one blue and 
meta l- poor \Fe/H] ^ —1.5 (jBrodie fc Strader 
20061 ). Peng et ^ ( 20061 ) found bimodal color dis¬ 
tributions for globular clusters in approximately 
100 early-type galaxies of the Virgo galaxy cluster. 
They found that, on average, red globular clusters 
account for 15% of the total population in faint 
{Mb ~ —15.7) galaxies and up to 60% in bright 
{Mb ~ —21.5 ) galaxies. Using the same obser - 
vational data. I^enda. Muriel fc Donzellil (200^ 
fitted power-laws to the projected number den¬ 
sity profiles, hnding typical logarithmic slope val¬ 


ues of a = —2.4 ± 0.2 and a = —1.7 ± 0.1 for 
red and blue populations, respectively. Although 
wit h a large scatter, other authors (see for exam¬ 
ple Lai:sen_et_aL 2001 : Bassino. Richtler fc Dirschl 


2008t iRichtler et al.N2012l : iD’Abrusco et al.N2014l) 
agree that blue GCs usually show a shallower dis¬ 
tribution than that corresponding to red GCs in 
the whole range of radii. Since the blue population 
is more extended than the red, it is expected that 
it will be more prone to be tidally disrupted than 
the red. Indeed, Sjv d ependence on clustercentric 
distan ce reported by iGoenda. Muriel fc Donzelli 
( 20091 ) has been detect ed only for blue but n ot for 
red globular clusters. iKartha et al. ( 20141) com¬ 
piled a sample of 40 galaxies and found that the 
relative fraction of blue to red globular clusters 
decreases with increasing galaxy environmental 
density for lenticular galaxies. They propose that 
the outer blue globular cluster population may be 
stripped away, giving a lower ratio of blue to red. 

From the theoretical point of view, an early 
work of Muzzio. Martinez fc Rabolli ( 19841) ana¬ 
lyzed the swapping of globular clusters in Virgo- 
like cluster of galaxies using ad-hoc N-body nu¬ 
merical simulations. These authors show that 
(smaller) galaxies in clusters can lose up to 30% of 
their globular cluster population, with half of this 
fraction ending up in the intracluster medium and 
the other half captured by other (larger) galaxies. 
Bekki fc Yahagil ( 2006h selected globular cluster 
particles in cosmological N-body numerical sim¬ 
ulations and found that intracluster globular clus¬ 
ters can contribute up to 40% of the total glob¬ 
ular cluster population in galaxy clu sters of mass 
1.0 to 6.5 X 10^‘^Mq. More recently. LSmith et al.l 
(|2013l) studied the impact of harassment on the 
dynamics of globular clusters in early-type dwarf 
galaxies of clusters. These authors modeled galax¬ 
ies that evolved for 2.5 Gyr within a potential field 
that mimics the effects of harassment in clusters, 
hnding that the dynamical behavior of the glob¬ 
ular clusters strongly depends on the fraction of 
bound dark matter remaining in the galaxy. They 
conclude that when 85% of the initial dark mat¬ 
ter halo mass has been removed, globular clusters 
begin to be stripped. 

Using N-body numerical simulations in the 
framework of the A Gold Dark Matter cosmo¬ 
logical model, we explore the tidal removal of 
globular clusters in a Virgo-like galaxy clusters, 
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focusing on the differences between red and blue 
populations. In Section [2] we present the numeri¬ 
cal simulation used in this work. In Section [3] we 
describe the method that we apply in order to se¬ 
lect dark matter particles from simulated galactic 
halos that mimic their globular cluster popula¬ 
tion. In Section |4] we present the main results 
about the removal of red and blue globular cluster 
populations. Finally, in Section [5] we discuss the 
main results and conclusions. 


2. The numerical simulation 


We address these issues using a dark matter- 
only cosmological numerical simulation that fol¬ 
lows the formation of a galaxy cluster-sized 
halo. This simulatio n was presented in detail by 
Ludlow et al. ( 201011 to which we refer the inter¬ 
ested reader for further technical details. In brief, 
a region destined to form a galaxy cluster at red- 
shift z = 0 is identified in a larger cosmological 
box of 100Mpc/i“^ comoving on a side. Then, 
this region is resimulated at a higher resolution 


using a zoom-in technique (jKlvDin et al.l 1200111 


while the remainder of the box is resampled at a 
lower resolution. The high-resolution region has 
1.44 X 10^ dark matter particles filling an amoeba¬ 
shaped volume of 1.39Mpc^h“^ at a starting red- 
shift of z = 19. The original simulation assumes 
a A Cold Dark Matter cosmology with the fol¬ 
lowing parameters: Hq = 73kms“^Mpc“^ (i.e. 
h = 0.73), Do = 0.25 and Da = 0.75, Both 
simul ations were performe d with the GADGET2 


code (ISuringel et al.N2005ri . High resolution dark 


matter particle masses of the resimulation are 
5.4 X 10^MQh~^ and the gravitational softening 
is 1.5Kpc/i“^ comoving. At redshift z = 0, the 
simulated galaxy cluster has a virial radius of 
Rvir = 1.1 Mpc h ~^, which corresponds to a virial 
mass Myir = 1.71 x IO^^'^Mq h~^. These values are 
compar able to those estimated for the Vi rgo Clus¬ 
ter (see Karachentsev fc Nasonov^ 2010l l. Within 


the virial radius (i.e. radius where the average 
mass density is 200 times the critical mass den¬ 
sity of the universe) at redshift z = 0 there are 
3.16 X 10® dark matter particles. This halo cor- 
responds to the hal o labeled ”hl4” in Table 1 of 
Ludlow et al. ( 2010ll which is isolated and virial- 
ized at z=0. In order to select the main cluster¬ 
sized halos and their g alactic subhalos, we use the 
SUBFIND algorithm ( Soringel et al. 2001) which 


identifies substructures inside friends of friends 
groups. We run SUBFIND on all 100 snapshots, 
equally spaced in log(a), from redshift z = 19 to 
z = 0 to follow the dynamical evolution of their 
globular cluster populations. In Fig. [1] we show 
the projected spatial distribution of dark matter 
particles (small gray dots) centered in the simu¬ 
lated galaxy cluster at redshift z = 0, where the 
solid circle shows the system virial radius. Using 
different colors we highlight dark matter particles 
selected as nine different globular cluster systems, 
following the method described in the next sec¬ 
tion. All the remaining globular cluster systems 
are depicted as gray dots. 


3. The model 


3.1. Globular cluster population 


Using SUBFIND (jSpringel et al.ll200ll) we iden¬ 
tified the 158 most massive dark matter halos that 
at redshift z = 0 are located inside the virial ra¬ 
dius of the galaxy cluster. At z=0, the less massive 
halo is resolved by at least 200 particles, which 
corresponds to M ~ 1 x 10 ^®Mq/i“^, while the 
most massive one has about 4000 particles, which 
corresponds to M ^ 2 x 10^^MQh~^. We follow 
the temporal evolution of each halo, computing its 
mass, density profile, orbit and accretion time i.e. 
the time when it enters the galaxy cluster virial ra¬ 
dius for the first time. In order to perform reliable 
statistics, we limited our sample to those halos 
that, at selection time, were resolved by at least 
700 particles. In order to maximize the amount of 
time that halos are orbiting the galaxy cluster, we 
exclude from our analysis those halos that have 
late accretion time (t > lOGyr). We have also 
excluded 19 halos that at the selection time (or 
even before) had a near neighbor halo. Our final 
sample has 38 halos. 


In Fig. m we show the clustercentric distance as 
a function of time for these 38 halos (gray lines), 
highlighting the nine halos also highlighted with 
the same color-coding in Fig. [TJ In order to 
minimize halos with possible tidal distorsions, we 
checked each halo during the infall process. We no¬ 
ticed that at accretion time many halos are already 
distorted by the clusters gravitational potential. 
However, two snapshots before that, halos show 
no evidence of such distortion. Gonsequently, we 
adopt this time as the snapshot to select globular 


3 































clusters. Selection time was, on average, 0.5Gyr 
before the accretion time (shown as vertical lines 
in Figl2]). We computed the mass density pro- 
file of halos at this time, fitti ng a NFW profile 
(jNavarro. Frenk fc White! 1996ll expressed as: 


PNFw{r) = 


Pnfw 


ir/rNFw)i^ + r/rMFwY 


( 1 ) 


where P%fw tnfw are the dark matter ha¬ 
los characteristic density and scale length, respec¬ 
tively. Typical accretion cosmic times range be¬ 
tween 4.0 and 9.3 Gyrs with a median of 5.2 Gyrs. 
Scale lengths range between 7.5 to 87.8Kpc/i“^ 
with a median of 10.5 Kpc h~^ while virial masses 
range between 5.3 x 10^° and 1.7 x IO^^Mq/i”^ 
with median 1.7 x 10 ^^Mq/i“^. In order to build 
up globular cluster populations of each halo, we se¬ 
lected a subsample of their d ark matter particles, 
using the m ethod outlined bviBullock fc Johnston 
(20^) and IPeharrubia. Navarro fc McGonnachie 
(|2008r) . In this method, a subset of particles is 
chosen to follow a given spatial mass density pro¬ 
file p(r) following its corresponding energy distri¬ 
bution function: 




(Pp dtp 


87r L Jq dtp^ Pe — tp 


\dtpJip=0- 


( 2 ) 


where tp is the relative gravit ational potential and 
£ is t he relative energy (see iBinnev fc Tremaine! 
198711 . 


In order to construct blue and red globular clus¬ 
ter populations for eac h halo at a ccreti on time 
we generate two analytic Hernauis^ (1993) density 
profiles 


PH{r) = 


Ph 


(r/r//)(l -h t/thY 


( 3 ) 


where is the scale length and p^ is the char¬ 
acteristic density. We scaled th = ^r^iFW adopt¬ 
ing 7 = 3 for blue and 7 = 0.5 for red globular 
cluster populations. With this choice for 7 at ac¬ 
cretion time, the projected density profiles at red- 
shift z = 0 are approximately power laws with 
slopes similar to those observed for blue and red 
populations (see Section [T]). At accretion tim e, we 
apply Equation (|2]) to the iHernauisd ( 199CJl and 
Navarro. Frenk fcWhitd (1 1996 1 profiles in order 
to construct their corresponding energy distribu¬ 
tions /ff(e) and fNFwY)-i respectively. 



Fig. 1.— Projected spatial distribution of dark 
matter particles (small gray dots) centered in the 
simulated galaxy cluster at redshift z = 0. The 
solid circle shows the system virial radius. Differ¬ 
ent colors show dark matter particles highlighted 
as blue globular clusters from nine halos, following 
the method described in Section |31 
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Fig. 2.— Glustercentric distances as a function of 
time for 38 halos (gray lines). Golors correspond 
to the nine halos shown in colors in Fig. [1] The 
black solid line is the virial radius of the cluster. 
Vertical lines show the time that each halo enters 
the cluster virial radius. 
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Then, we randomly select a fraction fn (e) / /at fw (e) 
of dark matter par ticles in bins of specific energy 
e + Ae. Following iBekki fc Yahasril ( 200ti[l we as¬ 
sume that the baryonic component is more concen¬ 
trated than the dark matter halo, so we truncate 
the selection at a cut off radius of rcuto ff= ^5o/3, 
where r^Q is the half-mass radius of the dark mat¬ 
ter halo. As an example, in Fig. [H we show the 
projected density profiles S(i?) of one selected halo 
as a function of the projected distance R, for the 
dark matter halo (black squares) and for the se¬ 
lected particles chosen to represent the blue (blue 
circles) and red (red triangles) globular cluster 
populations. In the top panel we show the profile 
at its accretion time (t = 4.9Gyrs,Zin = 1.25) 
and in the bottom panel at the present time 
{t = 13.7Gyrs,z = 0). This dark matter halo 
corresponds to that shown in light green in Fig. 

[U p and E The solid blac k line corresponds to 
a iNavarro. Frenk fc White ( 1996ll density pro¬ 
file with concentration c = 5.9 and virial mass 
Myir = 2.6 X Blue and red solid lines 

are not fits to the poin ts, but rather show the ana¬ 
lytic [Hernquis^ (1993) profiles chosen to represent 
both globular cluster populations. Dashed lines 
show power law fits to the points where we indi¬ 
cate the corresponding slope a for each globular 
cluster population at accretion time. We follow 
the temporal evolution of these blue and red pop¬ 
ulations in order to compare their properties with 
the corresponding populations observed in galaxy 
clusters and we present the main results in the 
next section. 


4. Results 

4.1. Density profile of globular clusters 

As can be seen in Fig. E the density profile of 
globular clusters is steeper at redshift zero than at 
the time of selection. We analyzed the evolution 
of these profiles as a function of redshift and found 
that the profiles become steeper immediately af¬ 
ter the halos enter the cluster. After that, slopes 
remain nearly constant independently of the evo¬ 
lution of the associated globular clusters. This 
indicates that, even in those cases where an im¬ 
portant fraction of globular clusters is removed, 
the remaining bound objects are redistributed in 
order to preserve a c ertain density profile. This 
is in agreement with iGoenda. Muriel fc Donzelii 



log(R) [Kpc ft-'] 


Fig. 3.— Projected radial density profile of one 
selected halo for dark matter particles (solid black 
squares), for blue (blue circles) and red (red trian¬ 
gles) globular clusters at redshift Zin {top panel) 
in which the halo enters the cluster’s virial radius, 
and at 2 = 0 {bottom panel). Top panel: The dark 
halo is modeled as a NFW profile (solid black line), 
while blue and red globular clusters are modeled as 
Hernquist profiles (solid blue and red lines, respec¬ 
tively). Vertical arrows show the Hernquist scale 
length rn for blue {ir^iFw) and red (0.5rAri?vv) 
globular cluster populations, and the NFW scale 
length {tfifw) for the dark matter halo. In both 
panels we show a power law fit to the data. Short 
dashed lines correspond to blue globular clusters 
and long dashed lines to red globular cluster pop¬ 
ulations. 
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( 2009f) . who found a lack of dependence between 
the slopes of the density profiles of globular clus¬ 
ters and the clustercen tric distance s in th e Virgo 
cluster. Nevertheless, iBekki et~n] ( 2003ll . using 
a simple numerical model, found that the radial 
number density of the globular cluster system be¬ 
comes steeper after the stripping of globular clus¬ 
ters. 


4.2. Globular cluster fraction evolution 


In order to follow the temporal evolution of 
globular clusters orbiting the galaxy cluster, for 
each halo, we estimate the fraction that is associ¬ 
ated with its own main halo at any given redshift. 
We assume that a globular cluster is associated 
with its main halo if its distance to the halo center 
is less than the most distant dark matter halo par - 
ticle identified by SUBFIND ( Sorineel et al. 2001 ) 
as part of this halo. In order to compute this frac¬ 
tion, at any given redshift z, we compute the ratio 
between the globular cluster number that we have 
selected at accretion redshift and the globular 
cluster number contained inside a sphere of ra¬ 
dius equal to the most distant dark matter parti¬ 
cle identified by SUBFIND. Although this selec¬ 
tion criterion may result in the inclusion of globu¬ 
lar clusters that are not dynamically linked to the 
halo, this criterion mimics what is usually done 
with the observational data, where globular clus¬ 
ters are selected up to a given radius. Moreover, 
we have found that our results are fairly robust if 
we adopt a binding energy instead of our distance 
criteria; average differences in the fractions plot¬ 
ted in figures 5, 6, 7 and 8 are always less than 
7%. 


Fig. [5] shows the fraction of globular clusters 
associated with its main halo as a function of time 
for both blue (top panel) and red [bottom panel) 
globular clusters (color coding for this figure is the 
same as that of Fig. [T]). The general trend is that 
at accretion time Un this fraction is, by definition 
f[zin) = 1 and decays to f[z = 0) ~ 0.6(0.75) at 
redshift z = 0 which means that about 40%(25%) 
of blue(red) globular clusters are removed from its 
halo to become part of the intracluster population. 
The steeper slopes in the upper panel of Fig. 
clearly show that blue globular clusters are more 
easily removed than red, as expected due to the 
shallower spatial distribution of the blue popula¬ 
tion (see Fig. SI). Some temporal oscillations are 
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Fig. 4.— Slopes of the projected density profiles 
of red and blue globular clusters at the redshift 
Zin in which they enter the cluster’s virial radius, 
versus the corresponding slopes at z = 0. Symbol 
and color types as in Fig. |3l Horizontal and long 
vertical arrows correspond to the median values of 
slopes for blue (solid blue line) and red (dashed red 
line) globular cluster populations. Short vertical 
arrows indicate observed slope values for red and 
blue CGs. Gray dashed line shows the identity 

C^20 ^Zin ■ 
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Fig. 5.— Associated fraction of globular clusters 
as a function of time for blue (top panel) and red 
(bottom panel) globular clusters for the same nine 
halos highlighted in Fig. [T] Color-coding as in 
FiglU 


noticeable in the fraction of globular clusters as¬ 
sociated with its halo (see for example the purple 
dashed line). Fig. [2]and[5]show that the local min¬ 
ima correspond to the successive passages of halos 
through the center of the galaxy cluster. However, 
it should be noticed that globular clusters are con¬ 
tinuously removed throughout the lifetime of halos 
not only due to pericentric passages. On the other 
hand, local maxima appear when their orbits in¬ 
tersect a stream of previously removed globular 
clusters. Also, we followed the fate of the stripped 
globular clusters and found that none of them are 
recaptured by other halos. 

4.3. Globular cluster fraction and orbital 
trajectory 

The strength of the tidal field generated by the 
gravitational potential of the galaxy cluster de¬ 
pends on the relative position of halos with re¬ 
spect to the cluster center. Thus, a correlation is 
expected between the fraction of globular cluster 
that remains associated with its halo and its or¬ 
bital parameters. Consequently, we searched for 
correlations between the fraction of globular clus¬ 
ters associated with its halo and the number of 
orbital revolutions completed by the halo from ac¬ 
cretion time tin to redshift z = 0. 

In Fig. 0 we plot the fraction of globular clus¬ 
ters at redshift z = 0 as a function of the number 
of orbits; points correspond to median values and 
vertical error bars are computed using the boot¬ 
strap resampling technique. There is a strong an¬ 
ticorrelation between the fraction of globular clus¬ 
ters at redshift z = 0 and the number of orbits. 
At least 30%(25%) of blue(red) globular clusters 
are stripped when halos have completed three or 
more orbits, while those that have completed only 
one orbit lose ~ 10% of their globular cluster pop¬ 
ulations. Based on their orbital trajectory, halos 
will cross the cluster center at different pericen¬ 
tric distances, experiencing tidal forces of differ¬ 
ent strengths that modify the fraction of globular 
clusters assigned. Fig. [7] shows the median frac¬ 
tion of globular clusters associated with each halo 
at z = 0 as a function of the minimun pericentric 
distance that each had throughout its whole tem¬ 
poral evolution. Vertical error bars are computed 
using the bootstrap resampling technique in bins 
of pericentric distance containing equal numbers of 
halos. Those that during their lifetime approached 
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the cluster core with pericentric distances less than 
150Kpc/i“^ lose approximately ^40%(25%) of 
their blue(red) globular cluster population, while 
halos that remain away from the core with peri¬ 
centric distances greater than 300 Kpc h~ ^ keep 
more than ~80%(90%) of their blue(red) popula¬ 
tion. Both Fig. El and [7] reinforce the interpre¬ 
tation that blue globular clusters are more easily 
removed from halos than red ones due to the dif¬ 
ferent slope of their density profiles. On average, 
at z = 0, halos have lost ^ 29% of their initial blue 
globular clusters; this number drops to ^ 16% for 
the red population, so it is expected that intra¬ 
cluster globular clusters will be mainly blue. 


The fraction of globular clusters that are associ¬ 
ated with its halo at redshift z = 0 shows a strong 
correlation with the number of orbits and the peri¬ 
centric distance; however, none of them are direct 
observables. In order to compare qualitatively our 
results with observations, we tested if the fraction 
of globular clusters associated, at any given red- 
shift, is a function of the clustercentric distance 
at that particular redshift. Fig. |8] shows the me¬ 
dian of the fraction of globular clusters associated 
with their halo as a function of the clustercentric 
distance for both red and blue populations. Dif¬ 
ferent line types correspond to different redshift 
ranges, averaged as indicated by the labels. We 
can see that associated clusters correlate with the 
clustercentric distance for three redshift ranges as 
shown. As expected, the effect depends on red¬ 
shift, being stronger at redshift zero. Although 
the fraction of globular clusters that are associ¬ 
ated with their halo is not a direct observable, 
th ese results can be compared with the findings 
of ICoenda. Muriel fc Donzellil ( 2009ll that the spe¬ 
cific frequency (Sn) oi the blue globular clusters in 
the Virgo Cluster increases with the clustercentric 
distance of galaxies. 


5. Discussion and conclusions 

Using a dark matter-only numerical simulation 
that follows the formation of a Virgo-like clus¬ 
ter of galaxies, we investigated the role of the 
cluster tidal field in the evolution of both red 
and blue globular cluster populations. At red¬ 
shift z = 0, the simulated cluster has a virial mass 
Myir = 1.71 X 10^^ Mq h~^, which is comparable 
to the mass of the Virgo Cluster. We selected 38 



Fig. 6.— Associated fraction of globular clusters 
as a function of the number of orbits at z = 0 
for blue and red globular clusters. Points cor¬ 
respond to median values and the vertical error 
bars are computed using the bootstrap resampling 
technique. Symbol and color types as in Fig. [H 



Minimum Ciustercentric distance [Mpc h'^] 

Fig. 7.— Associated fraction of blue globular clus¬ 
ters at z = 0 as a function of the minimum clus¬ 
tercentric distance. Points correspond to median 
values and vertical error bars are obtained using 
the bootstrap resampling technique. Symbol and 
color types as in Figure [3] 





























dark matter halos with M/ia/o > 1.08x 
that entered the cluster virial radius when the Uni¬ 
verse was younger than 10 Gyr. Using a Hern- 
squist profile, we randomly selected globular clus¬ 
ter particles in each dark matter halo, follow¬ 
i ng th e met hod outlined bv iBullock fc Johnston 
(20^) and IPeharrubia. Navarro fc McConnachie 
(|2008[ ). We used a Monte Carlo procedure to find 
the parameter that best reproduced the observed 
density profiles of both red and blue globular clus¬ 
ters. We summarize our principle results as fol¬ 
lows. 


- Blue globular clusters are lost more easily 
than red ones. This is an expected result 
that may be explained by the shallower spa¬ 
tial distribution of blue globular clusters. 
The stripping of globular clusters is a contin¬ 
uous process that increases with the succes¬ 
sive passages of halos through the center of 
the cluster. On average, halos have lost 16% 
and 29% of the red and blue globular clus¬ 
ter populations respectively, and therefore a 
significant number of globular clusters go to 
the intracluster medium. 


being stronger at redshift zero. This re¬ 
sult is in agreement with observational re¬ 
sults (see Coenda. Muriel fc Donz^ 20091: 


Forbes. Brodie fc Grillmaiij 19971) . 


• The density prohle slope of globular clus¬ 
ter systems becomes slightly steeper im¬ 
mediately after the halos fall into the galaxy 
cluster. Thereafter, it remains almost con¬ 
stant, independently of the evolution of 
the host halo. This result suggests that 
the density profile of globular clusters re¬ 
mains almost equal over the lifetime of ha- 
los , in agreement with the id e as su ggested 
bv ICoenda. Muriel fc Donzelli ( 2009| i. 


Our results confirm that the stripping of globu¬ 
lar clusters produced by the tidal field of clusters 
of galaxies is an efficient process that may result in 
significant loss of the original globular cluster pop¬ 
ulation. The magnitude of the effect will depend 
on the dynamical history of each galaxy. This phe¬ 
nomenon preferentially acts on the blue globular 
cluster population, and so it is expected that the 
intracluster globular clusters will be mainly blue. 


- None of the stripped globular clusters are 
recaptured by other halos. 

- Analyzing the bound fraction of globular 
clusters we found that at least 30% of blue 
globular clusters are stripped when halos 
have orbited the cluster of galaxies three or 
more times. For those halos that have com¬ 
pleted only one orbit the fraction removed is 
- 10%. 

- We found a strong anticorrelation between 
the associated globular clusters and the min¬ 
imum clustercentric distance that the ha¬ 
los have had throughout their history. Ha¬ 
los that have crossed the cluster core within 
150 Kpc from the cluster center during their 
lifetime, lost more than 40% of their blue 
globular cluster population, while halos that 
remain away from the core with pericentric 
distances greater than 300 Kpc, have associ¬ 
ated globular clusters greater than 80%. 
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this paper. We acknowledge Rory Smith, Ken 
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macros v5.2. 




Fig. 8.— Median of the fraction of the blue and 
red globular clusters associated with its parent 
halo as a function of the clustercentric distance at 
different redshifts: z = 0.06 (solid line), 2 : = 0.23 
(long dashed line) and 2 = 0.43 (short dashed 
lines). Color types as in Fig. [3] 
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